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© Method and apparatus for echo cancellation of a multicarrier signal. 

© When the transmitter and receiver of a duplex multicarrier modem are connected to a two-wire line, echoes 
of the transmitted signal are reflected back into the receiver by mismatches between the output impedance of 
the modem and the input impedance of the line. These echoes are cancelled by an echo emulator which 
operates on samples of the transmitted signal. The parameters of the emulator are adapted by correlating, in the 
frequency domain, the transmitted data with the received signal. A specialiyed phase-locked loop, which uses a 
unique phase detector, tracks both frequency offset and phase of the far echo seperately from the adaptation of 
the emulator parameters. Initial training of the canceller is performed with a simple cyclic training signal that 
allows the use of efficient fast transforms. 
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The present invention pertains to the transmission and reception of a multicarrier data signal, and more 
particularly, to the cancellation of the echoes of a transmitted signal that are superimposed on a received 
signal when both transmitter and receiver are connected to a bi-directional two-wire line. 

5 BACKGROUND OF THE INVENTION 

Duplex Transmission and Echo Removal: 

Many data communication systems require the simultaneous transmission and reception of data; this is 

70 called duplex transmission. In order to make the most efficient use of transmission equipment, it is often 
desirable that this duplex transmission be carried by a single bi- directional two-wire line. When, as is 
shown in Fig. 1, a transmitter 10 and a receiver 20 are connected to the same end of such a two-wire 
line-usually through a four-wire to two-wire connector or "hybrid" 30--a distorted version of the transmitted 
signal (an echo) is reflected back into the receiver. If the other end of the two-wire line is connected, 

/5 through a second hybrid 40, to a four-wire long-distance carrier system 50, then there may be another echo 
reflected by a third hybrid 60 at the far end of the four-wire line. These two echoes are called near-end and 
far-end echoes respectively; they would interfere with the detection of the signal from the far-end transmitter 
unless somehow removed. 

Early duplex modems such as those defined in the CCITT Recommendations V.22 and V.22 bis used 

20 separate frequency bands for the transmit and receive signals, and removed the echo(es) by filtering; this is 
called frequency division multiplexing, or FDM. More recent modems such as those defined in V.32 and 
V.32 bis, however, achieve more than double the data rate on a given line by using the full frequency band 
for the signals in both directions and removing the echo(es) by cancellation. 

The input impedances of the near-end and far-end two-wire lines to which a modem and a carrier 

25 system may be connected vary widely from line to line; furthermore, the impedance of a given line may 
also vary slightly with time. Therefore the characteristics of the echo(es) reflected back to the receiver vary 
widely, and the echo canceller must be adaptive; that is, upon initiation of a data connection it must learn 
the characteristics of the echo(es), and then track any small changes thereafter. 

30 Modulation Methods: 

The modulation method that has been most widely used for the transmission of high-speed data is 
Quadrature Amplitude Modulation or QAM (see, for example, J.A.C. Bingham, The Theory and Practice of 
Modem Design, John Wiley, 1988). This uses a single carrier at a frequency near the center of the usable 

35 frequency band; to distinguish it from the alternative method that will be described next, it will be referred to 
herein as single-carrier QAM (SCQAM). 

An alternative, and very efficient, method of transmitting high- speed data was described by Baran in 
U.S. Patent No. 4,438,511 "Packetized Ensemble Modem". The data to be transmitted are divided into 
blocks, and each block is further divided into sub-blocks of, typically, from one to eight bits. The bits in 

40 each sub-block are used to modulate one of a set of carrier signals, which are evenly spaced across the 
usable frequency band. All the modulated carriers are then added together, and the composite signal is 
transmitted as a "multicarrier" signal. 

PRIOR ART 

45 

Echo Cancellers: 

Echo cancellers (ECs) for SCQAM may be either signal-driven or data-driven, and may operate at either 
the symbol rate of the transmitted data or the Nyquist rate of the channel. All combinations of these options 

so have been described in the literature, and several patents have been granted (see, for example, S.B. 
Weinstein, U.S. Patent No. 4,131,767 and K.Mueller, U.S. Patent No. 4,087,654). 

Most of the early ECs were signal-driven as shown in Figure 2(a); the transmit signal from the 
Encoder/Modulator 70 is input to the hybrid 30 and aiso-usually in the form of digitized samples at the 
Nyquist rate-to an Echo Emujator 80 in the form of a Finite Impulse Response (FIR) filter (also known as a 

55 Tapped Delay Line, or TDL). The output of the emulator is then subtracted 90 from the output of the hybrid 
(comprising both the far-end signal and the echo of the local transmitter), and the difference is passed to 
the receiver. This difference is also correlated with the samples of the transmit signal stored in the emulator, 
and the products are used to update the parameters (tap values) of the emulator; this is known as Least 
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Mean Squares (LMS) adaptation, and it should converge to a state where the echo component of the 
received signal is almost completely cancelled. The main disadvantage of this method is that the spectral 
characteristics of the (necessarily) over-sampled transmit signal are such that adaptation of the tap values is 
very slow, and poorly conditioned numerically; this is known as the eigenvalue-spread problem, and it is 

5 often encountered in LMS adaptation. 

More recent ECs for SCQAM have been data-driven as shown in Fig. 2(b); the baseband data is 
encoded 100, and input to the Echo Emulator 110, which emulates both modulator and echo path. Because 
the data are scrambled (i.e., uncorrelated with each other) there is no spread of eigenvalues, and adaptation 
is much better behaved. Furthermore, in order to simplify the synchronization of the receiver. with the clock 

w in the far transmitter, both the data and the output of the hybrid are sampled at the Nyquist rate of the 
channel. 

Training: 

15 Unless some special measures are taken (see, for example, D.D.Falconer, "Adaptive Reference Echo 
Cancellation", IEEE Trans. Commun., pp. 2083-2094, Sept. 1982), the signal from the far-end transmitter 
will act like noise as far as the training of the EC is concerned, and will seriously degrade and/or slow down 
the adaptation. Therefore, most ECs are trained in two stages. During the first, half- duplex, stage the far- 
end transmitter is silent, and the received signal comprises only the echo(es); in the second, full-duplex, 

20 stage the received signal comprises both echo(es) and far-end signal, and the step-size of the adaptation 
must be greatly reduced. 

Training during the first "stage is often done using a repetitive (i.e., cyclic) signal such as a chirp 
sequence (see D.C. Chu, "Polyphase Codes with Good Periodic Correlation Properties", IEEE Trans. Inf. 
Theory, pp. 531-532, July 1972). The echo is transformed into the frequency domain and its transform is 

25 divided by the known transform of the transmit signal. Ling and Long ("Correlation-Based Fast Training of 
Data-Driven Nyquist In-Band Echo Cancellers", IEEE Intl. Conf. on Commun., pp. 1280-1284, April 1990) 
proposed that the period of the sequence should depend upon the delay of the far echo; such a period 
would not, in general, consist of a number of samples that was an integer power of two, and therefore the 
transformations could not be performed by the efficient Fast Fourier Transform; the much less efficient 

30 Discrete Fourier Transforms (DFT) would have to be used. 

Near and Far ECOs: 

If the round-trip propagation delay through the four-wire part of a connection is longer than the duration 
35 of the near-end echo, then the near and far echoes will be separated by a quiet period. Most ECs are 
therefore implemented as two separate cancellers (NEC and FEC) separated by an adjustable bulk delay. 

Because the far echo traverses one or more carrier systems it may have more impairments than the 
near echo; in particuar. it may be shifted slightly in frequency. The FEC must therefore be able to learn and 
track a small frequency offset between the transmitted signal and its far echo. Methods of doing this have 
40 been described by several authors (e.g., R.D. Gitlin and J.S. Thomson, "A Phase Adaptive Structure for 
Echo Cancellation", IEEE Trans. Commn., pp. 1211-1220, Aug. 1978). 

Multicarrier Modulation: 

46 Multicarrier modulation systems can be divided into two categories: (a) those that use filtering to 
achieve a degree of orthogonality between the signals in the sub-channels, and (b) those that modulate the 
carriers with unfiltered signals. Since the carriers in the first category are modulated and those in the 
second are keyed, these systems have been designated as MCQAM and MCQASK, respectively. The latter 
require much less computation because they can use very efficient Fast Fourier Transforms for modulation 

so and demodulation; furthermore they can use many more (mere closely spaced) carriers, and thereby 
achieve greater immunity to both channel distortion and impulse noise. 

Echo cancellation of MCQAM signals has been discussed by N.H. Lewis et al in "An Efficient Data 
Driven Echo Canceller for Multicarrier Modulation Systems" (IEEE Trans. Commun., in review), but neither 
the amount of echo suppression nor the speed of adaptation reported were adequate for high-speed data 

55 transmission over a wide range of channels. 

Echo cancellation of MCQASK signals has been discussed briefly by J.A.C. Bingham in "Multicarrier 
Modulation for Data Transmission: an Idea Whose Time Has Come" (IEEE Communications Magazine, May 
1990). As pointed out there, the distortion in the echo path will cause the data modulated onto each carrier 
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to produce echoes in many sub-bands besides its own; there will be both intra-channel and inter-channel 
4 echoes. Consequently, a straightforward implementation of a data-driven EC would require a prohibitive 

amount of computation. On the other hand, a conventional signal-driven canceller that operated on the 
composite multicarrier signal would suffer from the eigenvalue-spread problem. 

5 

SUMMARY OF THE INVENTION: 

In order to overcome the aforementioned problems the invention provides an improved method and 
apparatus for echo cancellation of a multicarrier signal which method is evident from independent claims 
w 1,8,9,10, and 11 and which apparatus is evident from independent claims 12,18,19,20, and 21. Further 
advantageous features, aspects, and details of the invention are evident from the dependent claims, the 
description and the drawings. The claims are intended to be understood as a first non-limiting approach of 
defining the invention in general terms. 

The full process of echo cancellation can be conveniently considered as three separate operations: the 
75 emulation of the echo, subtraction from the output of the hybrid, and adaptation of the parameters of the 
emulator. These operations can be performed in the time domain or the frequency domain, and the 
emulation and adaptation can be performed on either the transmit data or the transmit signal. 

One aspect of the present invention is to resolve the signal-driven/data-driven dilemna and thus avoid 
the problems associated with each approach. This is achieved by choosing the most appropriate domain 
20 and operand for each of the three operations. 

The emulation is performed upon the transmit signal; this avoids the problems of multi-dimensional 
emulators that would be needed for multicarrier signals if the emulators were data driven. The emulation 
can be performed in either the time domain using convolution or in the frequency domain using overlap- 
save transforms, both techniques are well-known in the art. 
25 The subtraction of the emulated echo must, however, be performed in the time domain. This preserves 
independence of the timing in transmitter and receiver, and also avoids the problems of the non-cyclic 
nature of the transmit signal. 

The adaptation of the parameters of the emulator is performed in the frequency domain by correlation 
with the transmit data. This avoids the problems of eigenvalue spread, and makes the adaptation process 
30 very well conditioned numerically. 

A second aspect is to do as much as possible of the signal processing in a block mode in the 
frequency domain (rather than sample by sample in the time domain) in order to take advantage of the 
computational efficiencies of the Fast Fourier Transform (FFT). 

35 BRIEF DESCRIPTIONS OF THE DRAWINGS 

Figure 1 shows the prior art connection of a transmitter and a receiver of a modem to a two-wire line, 
which is connected at its other end to a four-wire carrier system. This, in turn, is teminated by another two- 
wire line, which is terminated at its far end by the input impedance of the far modem. 
40 Figure 2a shows a signal-driven echo canceller of the prior art. 

Figure 2b shows a data-driven echo canceller of the prior art. 

Figure 3 is a detailed depiction of one embodiment of the present invention for the cancellation of a 
near-end echo. 

Figure 4 shows an alternative, prior art, embodiment of an echo emulator. 
45 Figure 5 shows the preferred method of compensating for frequency offset in the far echo. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Before beginning to train the echo cancellers the delay, T d) and the frequency-offset, f 0( f, of the far echo 
so must be measured. Methods of measuring the delay and the offset have been described in CCITT 
Recommendation V.32 and by J. Cioffi ("A Fast Echo Canceller Initialization Method for the CCITT V.32 
Modem", IEEE Trans. Commun., pp.629-638 t May 1990), respectively. 

After the delay has been measured and an initial estimate of any frequency-offset has been made there 
are three modes of operation for the canceller of the present invention: 
55 (1) initial frequency-domain training using cyclic signals, 

(2) adaptation of the canceller parameters in a half-duplex mode (that is, without double-talk), and 

(3) continual slow adaptation of the parameters during double-talk. 

In the first embodiment of the invention the only difference between operation in modes (2) and (3) is in 
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the adaptation step size; during double-talk it must be much smaller in order to filter out the part of the 
received signal that comes from the far-end transmitter, which, paradoxically, acts like noise for the echo 
canceller. Operation in modes (1) and (2) will be described in reverse order, as is common practice in the 
art. 

5 

Generation of the Transmit Signal: 

One method of calculating a block of samples of an MCQASK signal and transmitting them sequentially 
onto a line has been described by Keasler et al in U.S. Patent No.4,206,320 and by Baran in U.S. Patent 
10 No. 4,438,51 1 . The steps are shown in blocks 200-230 of Fig. 3, and can be summarized as follows: 

(1) the serial input data are assembled into blocks using a serial-to-parallel converter 200, and each 
block is divided into N c sub-blocks as described by Hughes-Hartogs in U.S. Patent No. 4,679,227, 

(2) each block and constituent sub-blocks are encoded 210, and the binary data in each of the sub- 
blocks are converted into two-dimensional (i.e., complex), multi-valued signals, which are to be used to 

is modulate the carriers, 

(3) each set of N c complex, multi-valued signals (ignoring, for the moment, the Augmentation operation 
215) is then used to modulate N c carriers at frequencies evenly spaced from ni f to n 2 f, where N c = r\ 2 
- ni +1. This modulation is effected by an Inverse FFT (IFFT) 220, which converts the frequency- 
domain samples into a set of N real, time-domain samples. In the first implementation Af was 7.8125 

20 ( = 8000/1024) Hz, and N was 1024; N c varies with the transmission channel used, but is typically about 
380. 

(4) the last L samples of each block are also prefixed 230 to the front of that block in order to eliminate 
inter-symbol interference in the far-end receiver by providing a guard period during which the propaga- 
tion transients of the line may decay, and 

25 (5) the resultant (N + L) samples are converted, in parallel-to- serial converter 240, into a sequence of 
signals 250, which are passed to a digital-to-analog converter 260, a low-pass filter 270, a hybrid four- 
wire to two-wire connector 280, and thence to the two-wire line. 

Augmentation of the Transmit Signal for Improved Correlation: 

30 

For the most precise adaptation of the canceller it is desirable that the set of frequency-domain samples 
that is input to the correlator 290 and the IFFT 220 have a slightly wider bandwidth than just the N c carriers 
that are used for transmission of data. Therefore the set of N c complex, multi-valued samples may be 
augmented, in 290, by a few fixed (i.e., not changing from block to block) samples, which will result in 
35 unmodulated carriers below ni f and above r\ 2 f. The total number of carriers used will be defined as N t01 . 

Emulation and Subtraction of the Echo: 

Fig. 5 shows the full implementation of the invention, with cancellation of both near and far echoes. For 
40 the sake of clarity in a first explanation, however, Fig. 3 shows only the cancellation of a near echo. 

The signal received from the two-wire line, which is output from the hybrid 280, is passed to an anti- 
aliasing low-pass filter 340, and an analog-to-digital converter, ADC 350. The digital samples of this signal 
390 are then input to the plus terminal of subtractor 380. During duplex data transmission and reception 
("doubletalk") this signal is the sum of the signal from the far transmitter, the near echo, and, sometimes, a 
45 far echo. In order, however, to simplify the first explanation of the echo canceller it will be assumed that the 
far transmitter is silent, and that also there is no far echo. 

The sequence of transmit signals 250 is input to near echo emulator 320, and the output of this 
emulator is applied to the minus input of subtractor 380. Thus, a near echo is cancelled by subtracting from 
the received signal a locally generated-that is, emulated-echo. 
so En one embodiment of the invention, the echo emulator 320 is implemented as a tapped delay line, 
which performs a conventional convolution in the time domain: 
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where the u k and v k are the samples of the input to and output from the emulator, and the h m are the most 
recent estimates of the samples of the truncated impulse response of the echo path, which are derived from 
the output of the IFFT 300. 

In an alternative embodiment of the present invention echo emulation may be performed in the 
5 frequency domain as shown in Fig. 4, where the Echo Emulation block 320 has been expanded to show the 
following operations. 

(a) Samples of the transmit signal 250, defined as Xj, are grouped into blocks j = kN + 1 to (k + 1)N, and 
transformed by FFT 400. They are also delayed by half a block (i.e., by N/2 samples) in 410, and 
samples j = (k + 1/2)N + 1 to (k + 3/2)N are transformed by FFT 420. 
;o <b) The two frequency-domain transforms of the transmit signal are each multiplied in 430 and 440, sub- 
channel by sub-channel, by the frequency-domain parameters of the echo, 

(c) The two sets of products are transformed back into the time domain by IFFTs 450 and 460, and 

(d) Because of the non-cyclic nature of the transmit samples and the duration of the echo impulse 
response, each of the sequences of time samples output by 450 and 460 is valid during only part of the 

75 time. Therefore switch 470, in the position shown, selects the output of 450 for j = (k + 1/4)N + 1 to 
(k + 3/4)N; it then switches to the other position and selects the output of 460 for j = (k-1/4)N + 1 to 
(k + 1/4)N. 

Adaptation of the Parameters of the Echo Emulator. 

20 

The samples of the modified received signal (received signal minus emulated echo) are grouped into 
blocks by the Serial-to- Parallel converter 370 and then input to the FFT operation 360, and transformed into 
N t0! complex, frequency-domain samples, Y\. These and the frequency-domain samples of the transmit 
signal Xn, are input to an LMS algorithm, comprising Subtraction 330, Correlation and Updating 290, IFFT 
25 300 and FFT 310, which stores and updates the estimates H n (n = 1 to N tot ) of the frequency-domain 
parameters of the echo path. A set of N tot error samples, E 0( is~calculated by 

E n = Y; n -H n X n (2) 

30 and a set of temporary, updated estimates of the H n are calcuated according to 

H*„ = H n -5E n X* n (3) 

where the * denotes the complex conjugate, and 8 is a small step size chosen, as is well known in the art, 
35 for the best compromise between noise averaging and speed of adaptation of the estimates. 

These temporary estimates H' n are then transformed into the time domain by the IFFT 300, and all but 
the first Mi time samples are set to zero. The new, truncated set of h m for m = 1 to Mi is then input to the 
echo emulator for use with the next block of transmit signals. The truncated set is also retransformed back 
into the frequency domain by FFT 310 in order to generate a new, smoothed set of estimates H n . 

40 

Cancellation of Far-End Echoes: 

The method of cancelling far echoes is the same as that for near echoes except for the addition of three 
bulk delays-two serial and one block-and an algorithm to compensate for frequency offset; these features 

45 are shown in Fig. 5, where several of the individual operation blocks of Fig. 3 have been combined into 
composite blocks (e.g., 290, 300, 310 and 330 are shown as F-D SC&U 530). 

Because of the propagation delay in the four-wire path (the main component of the round-trip delay T d ) 
shown in Fig. 1, there may be a significant period of time between the transmission of a signal and the 
return of a far echo to the receiver. It is undesirable (in terms of both the amount of computation required 

so and the accuracy that is achievable) to use enough parameters in the far echo emulator to span this quiet 
period, so nearly all prior art echo cancellers for SCQAM use a bulk delay, during which the data and/or the 
transmit signal is stored but is not expected to produce an echo. For the MCQAM echo canceller of the 
present invention this function is performed by the Serial Bulk Delay 520, which delays the transmit signal 
byT 8l . 

55 The block of encoded, complex data samples plus a few constant samples, which together are to be 
used to modulate the carriers, are input, as described previously, to the IFFT and Cyclic Prefix 500 and the 
F-D SC&U 530 and now also to the Block Bulk Delay 510; the output of this delay is input to a second F-D 
SC&U 540. Similarly, the samples of the transmit signal, which are output from the Parallel to Serial 
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converter (P/S) 240, are input to Serial Bulk Delay 520, and the output of this delay is input to the Far Echo 
Emulator 550, which, in the first embodiment of the invention, was implemented in the time domain, but 
might, as already discussed, alternatively be implemented in the frequency domain using overlap-save 
FFTs and IFFTs. The output of the emulator is then input to Phase Splitter 560. 
5 The Phase Splitter 560 and Multiplier 600 rotate the emulated echo by a continually changing phase 
angle if the echo has been subjected to a frequency offset. The output of the multiplier is subtracted from 
the modified received signal in Subtracter 610. 

The total delay of the transmitted signal through the serial bulk delay (T s1 ) and phase splitter (T ps ) must 
be equal to T d , the round-trip delay of the echo, so ideally 

T.i = T d - T ps (4) 

[In practice, T s1 is set slightly less than this in order to allow for any error in measuring T d .] 

The parameters of the Far Echo Emulator are then adapted in the same way as those of the Near Echo 
75 Emulator, using the second serial-to- parallel converter S/P 660, FFT 650, and F-D SC&U 540. 

In order, however, that the same block timing can be used for the input to the far echo F-D SC&U as for 

the transmit signal and the near echo F-D SC&U, it is necessary to delay by T^, in Serial Bulk Delay 670, 

the input to the second serial-to parallel converter S/P 660 so that the total delay of both echo and emulated 

echo is an integer number of block periods. That is, T^ is defined by 

20 

T d + Tsa = K(N + L)T (5) 

Then, in order that blocks of the transmit signal be correlated with the appropriate blocks of the echoes, 
the Block Bulk Delay must be equal to K blocks (i.e., T b = K(N + L)T). 

25 

Adaptation for Frequency Offset Compensation: 

The F-D SC&U 540 continuously learns the complex frequency-domain parameters of the far echo, Hf„ 
which it passes to the Far-Echo Emulator 550, and the algorithm dO 570. This last unit calculates the 
30 change of the phase of each parameter and the average change over ail N c parameters. 

By comparing the smoothed parameters, Hf in , calculated after processing of the ith block with the 
previous set Hf M , n , set of smoothed changes can be defined by 

dHf iin = HfLn-Hfi-Ln (6) 

35 

Then the small angle by which each complex parameter has rotated is given aproximately by 

d^,, n = lm(^i, n .d^i, n )/lH_f iin .Hf, 1( n: 2 (7) 

40 where the * denotes the complex conjugate, and the small angle approximation that dO/ a arctan[d$] has 
been used. 

The division in equation (7) is a tedious operation in most digital signal processors, so it is dispensed 
with. The result is that the contributions from each frequency sub-channel to the estimate of 60/ are 
weighted by the power received in that channel--a not unreasonable approximation. 
45 An average phase change over all N c sub-channels is then calculated by 

d^ « ~ E d0. n (8) 

so N c n^n, 

This average phase change is input to a Low-Pass Filter 580. A filter for a conventional Phase-Locked 
Loop (PLL), would, as is well understood in the art, have a transfer function 

55 
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9. S z[z - (l-B) ] 

F(2) = = (9) 

5 d0. (z - l) 2 

where 6 and £ are small quantities chosen for the best compromise between transient response and noise 
filtering. The input to this filter is however, a phase change, so the output is also; to calculate the average 

10 phase of the echo path another integration 620 must be added. 

The angle 6\ that is calculated in this way is the average, across all frequency bands and across the 
period of a block, of the phases of all of the multicarrier echo path sub-channels; if there is frequency offset 
the actual phase of each of the sub-channels will, of course, change linearly with time across the period. 
Therefore the phase that is input to the rotator 600 for sample k (k - 1 to N + L) should have two 

75 components: 

6 ]M = e ( + 2TTkf oH(i (10) 

where f offii is the estimate, after the ith block, of the frequency offset. This estimate is proportional to the 
20 filtered phase change, so equation (10) can be rewritten as 

hk = *i + W + L)](fl,.y (10a) 

and implemented by multiplier 630 (with 8 ~ k/(N + L)) and summer 640. 
25 The angle Q\ is then used in 590 to generate the complex number e* 0 = cos B + jsin d. A look-up table 
was used in the first embodiment, but various methods of calculation that are well-known in the art might 
alternatively be used. The complex number is then input to the multiplier 600 which generates just the real 
part of the emulated and rotated far echo. 

30 initial Frequency-Domain Training 

Ling and Long (loc. cit.) showed that if the round-trip delay T d is such that 
(K'N-M 2 )T< = T d < (K'N + M,)T (11) 

35 

where K' is any integer, then, because of the cyclic nature of the training signal (period = NT), near and far 

echoes will appear to overlap, and will not be completely separable. 

Their solution to this problem was to adust the period of the training signal according to T d . The 

algorithm for doing this was fairly complicated, and, because the resultant period did not necessarily contain 
40 a number of samples that was a power of two, the computationally-efficient fast algorithms could not be 

used for the transformations. 

The solution to this problem that is incorporated into the present invention is to use a period for the 

training signal which (a) is equal to the block length used during data transmission, (b) comprises a number 

of samples that is an integer power of two. (c) is independent of the delay T d , and (d) is longer than the 
45 longest T d encountered on terrestrial transmission links. 

Calculation of Echo Path Impulse Responses: 

The preferred cyclic transmit signal is a chirp sequence with a fixed period of 128 ms. The 512 
so complex samples (at 4000 samples/s) are modulated in sequence onto all 512 carriers at frequencies from 0 
to 4.0 kHz using the IFFT 220, and transmitted as described previously. The method of calculating the echo 
path response is slightly different for terrestrial and satellite links. 

Terrestrial Links: 

55 

For connections that do not involve transmission via a satellite the delay of the last part of the far echo, 
(T d + M 2 T), is less than the symbol period, NT (= 128 ms), of the preferred embodiment. Therefore the 
echo path response is calculated as follows. 
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(a) The transform of the received signal is divided by that of the transmitted chirp sequence. Since the 
samples of the chirp sequence are all of equal magnitude, which can be normalized to unity, division can 
be effected by multiplying by the complex conjugate of the transmitted sample. 

(b) The dividend is IFFTed into the time domain, and the samples separated into two sets-those in the 
5 ranges of t = 0 to Mi T and from T d to T d + M 2 T--and the rest are set to zero. 

(c) The near echo response is averaged over several blocks, and loaded, as an initial setting for 
subsequent adaptation, into the F-D SC&U 530. 

(d) The far echo impulse response is derotated in accordance with the learned frequency offset, and a 
best estimate of the frequency- domain and time-domain parameters of the far echo is made after each 

70 block. 

(e) The far echo response, calculated in this way, is averaged over several blocks, its delay is reduced 
by T d »so that the echo appears to start at the beginning of a block-and the resultant frequency-domain 
and time-domain parameters are loaded into F-D SC&U 540 and Far Echo Emulator 550. 

(f) The parameter, K, of the Block Bulk Delay is set, as indicated by equation (5), to one. 
75 (g) The estimate of frequency offset is updated with each block after the first. 

Satellite Connections: 

For these the round-trip delay is greater than 450 ms. Therefore at least three cycles of the chirp 
20 sequence are transmitted before the far echo arrives at the receiver; during that time the received signal 
390 comprises only the near echo. The near echo path response can therefore be almost completely 
learned before the far echo arrives. 

For the Kth cycle (the cycle during which the first part of the far echo will arrive) and all subsequent 
cycles of the initial training signal, the near and far echoes may appear to overlap, and the estimate of the 
25 near echo cannot be further improved. The parameters are therefore fixed at those learned from the first (k- 
1) cycles, and the emulated near echo subtracted from the received signal. This emulation and subtraction 
is performed in the time domain in the first embodiment using 320 and 380, but it might alternatively be 
performed in the frequency domain. In either case, the time-domain difference signal, which is assumed to 
comprise only the far echo, is input to the far echo initialization algorithms, which are thereafter the same as 
30 for terrestrial connections. 

Claims 

1. A method of cancelling the echoes of a signal transmitted via a two-wire line in a data modem 
35 comprising a multicarrier transmitter which transforms blocks of encoded frequency-domain samples 

into blocks of time-domain samples and a multicarrier receiver which performs the inverse operation, 
said transmitter and receiver being both connected to said two-wire line, which echoes are reflected 
back into said receiver by an imperfect match between the output impedance of said modem and the 
input impedance of said two-wire line, said method comprising the steps of: 
40 (a) forming a first estimate of the frequency-domain characteristics of the path traversed by said 

echo; 

(b) transforming said estimated frequency-domain characteristics of step (a) into the time-domain to 
form an estimate of the impulse response parameters of said echo path; 

(c) calculating the expected time-domain samples of the echo of said transmitted signal based upon 
45 said estimated impulse response parameters of step (b); 

(d) subtracting said expected echo signal of step (c) from the actual signal received from said two- 
wire line in order to calculate a modified received signal; 

(e) transforming a first block of time-domain samples of said modified received signal of step (d) into 
the frequency domain; 

so (f) correlating said transformed modified received signal of step (e) with said encoded frequency- 

domain samples of said transmitted signal; 

(g) calculating a new set of estimated frequency-domain echo path characteristics by decreasing 
each of said first estimated frequency-domain characteristics by a small quantity that is proportional 
to said correlation of step (f); and 
55 (h) mproving the estimate of said frequency-domain characteristics by repeating steps (a) - (g) using 

subsequent samples of said transmitted and received signals. 

2. The method of claim 1 in which said transformations from time domain to frequency domain are 
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performed using a discrete Fourier transform, and said transformations from frequency domain into 
time domain are performed using an inverse discrete Fourier transform. 



The method of claim 1 in which said transformations from time domain to frequency domain are 
performed using a fast Fourier transform, and said transformations from frequency domain to time 
domain are performed using an inverse fast Fourier transform. 

The method of one of the preceding claims in which said estimates of the characteristics of the echo 
path comprise just one complex number for each of the carriers used in said transmitted signal. 

The method of one of the preceding claims in which said encoded frequency-domain samples occupy 
a first frequency band, and are augmented by some other samples that are not used to convey data 
information, so that the composite set, which is transformed into said set of time-domain samples and 
is also correlated with said modified received signal, occupies a second frequency band that is wider 
than said first frequency band. 

The method of claim 1 in which the number of parameters of said estimate of the impulse response of 
the echo path is less than the number of parameters of said estimate of the frequency-domain 
characteristics, and a smoothing of said frequency-domain parameters is effected by the steps of: 
(i) changing said frequency-domain characteristics in order to generate a first set of new characteris- 
tics; 

(j) transforming said first set of new characteristics into the time-domain, and truncating the time- 
domain impulse response in order to generate a new set of parameters of said estimated impulse 
response; and 

(k) transforming said new, truncated estimated impulse response back into the frequency domain to 
generate a second set of new frequency-domain characteristics. 

The method of claim 6, in which said transformations of steps (j) and (k) are performed using Inverse 
Fast Fourier Transforms and Fast Fourier Transforms, respectively. 

A method of cancelling a far echo of a transmitted signal in a data modem comprising a multicarrier 
transmitter which transforms blocks of encoded frequency-domain samples into blocks of time-domain 
samples and a multicarrier receiver which performs the inverse operation, said transmitter and receiver 
being both connected to a communication link that comprises, in turn, a first two-wire line, one or more 
four-wire terrestrial carrier systems, and a second two-wire line which echo is reflected back to said 
receiver from said second two-wire line, and which may be offset in frequency, said method comprising 
the steps of: 

(a) forming first estimates of the round-trip delay of said far echo, and of the frequency-domain 
characteristics and frequency offset of the path traversed by said far echo; 

(b) transforming said estimated frequency-domain characteristics of step (a) into the time-domain to 
form an estimate of the impulse response parameters of said far echo path; 

(c) calculating the expected time-domain samples of the echo of said transmitted signal based upon 
said estimated impulse response parameters of step (b); 

(d) rotating said expected time-domain samples of step (c) by an angle that increases linearly with 
time in proportion to said estimate of the frequency offset; 

(e) subtracting said rotated expected echo signal of step (d) from the actual signal received from 
said two-wire line in order to calculate a modified received signal; 

(f) delaying said samples of said modified received signal of step (e) by a time that is equal to the 
difference between an integer number of periods of said blocks and said round-trip delay of step (a); 

(g) transforming a first block of time-domain samples of said modified and delayed received signal 
of step (0 into the frequency domain; 

(h) delaying said blocks of said encoded frequency-domain samples by a time that is an integer 
multiple of the period of said blocks; 

(i) correlating said transformed modified received signal of step (g) with said encoded and delayed 
frequency-domain samples of step (h); 

fl) calculatng a new set of estimated frequency-domain echo path characteristics by decreasing each 

of said characteristics by a small quantity that is proportional to said correlation of step (i); 

(k) calculating a set of estimates of the phase angles by which each of said new frequency 
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Direct-Vent Safety 

For clean, comfortable heat 
when and where you want it. 
The gas flame is completely 
sealed inside, away from all 
contact with room air. Air 
necessary for combustion is 
drawn into the chamber from 
outside through direct wall 
vent, and all combustion 
products are vented directly 
to the outside. Only clean, 
gentle heat, in continuous 
circulation is delivered into 
the room. Windows and 
doors can stay closed 
because no room air is used 
for combustion. 



"Matchless" Pilot Ignitor 

All models include a Piezo pilot ignition and a glass observation door. 
Assures convenient pilot lighting without matches. 

Automatic Thermostats 

(No Electricity Required — Included with all Models) 
Available with your choice of a self-generating millivolt wall thermostat or 
a self-contained bulb thermostat with a full range of adjustment on the 
gas valve dial. All models feature 100% safety shut-off as standard 
equipment. Mobile home models equipped with wall thermostat only. 

Contemporary Finish 

Light, powder-coat scratch-resistant finish to blend with any decor. 

Conversion Kits Available 
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characteristics has rotated; 

(I) averaging said changes of phase angles of step (k) over the set of said characteristics; 
(m) filtering said average phase angle of step (I) to generate an improved estimate of the frequency 
offset of said far echo; and 
5 (n) repeating steps (b) through (m) in order to continually improve said estimates of step (a). 

9. A method of initially estimating the parameters of near and far echoes of a transmitted signal in a data 
modem comprising a multicarrier transmitter which transforms blcoks of encoded frequency-domain 
samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 

w operation, said transmitter and receiver being both connected to a terrestrial communication link that 
comprises, in turn, said first two-wire line, one or more four-wire terrestrial carrier systems, and said 
second two-wire line, which echoes are reflected back into said receiver by said first and second two- 
wire lines respectively, said method comprising the steps of: 

(a) repeatedly transmitting a cyclic signal that has a period that is equal to the duration of said sets 
is of encoded frequency-domain samples, and longer than the delay of said far echo; 

(b) transforming into the frequency domain sets of samples of the signals received via said echo 
paths; 

(c) dividing said transform of said sets of samples of step (b) by the transform of said cyclic signal 
of step (a) to form a first estimate of the combined frequency-domain characteristics of both of said 

20 near and far echo paths; 

(d) calculating the time-domain parameters of said near echo path by transforming said first estimate 
of step (c) into the time domain, and setting to zero all samples outside the span of said near echo; 

(e) calculating a first set of time-domain parameters of said far echo by setting to zero all samples of 
said transformed first estimate of step (d) which are outside the span of said far echo; 

25 (f) generating a second set of time-domain parameters of said far echo by rotating said first set of 

time-domain parameters of step (e) by the angles needed to compensate for the estimated 
frequency offset of said far echo path; 

(g) transforming said second set of time-domain parameters of step (f) into the frequency domain; 
and 

30 (h) for all sets of samples of said actual received signal after the first, averaging the change of phase 

angle of said frequency-domain parameters and thereby determining an update of the estimated 
frequency offset. 

10. A method of initially estimating the parameters of near and far echoes of a transmitted signal in a data 
35 modem comprising a multicarrier transmitter which transforms blocks of encoded frequency-domain 

samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 
operation, said transmitter and receiver being both connected to a communication link that comprises, 
in turn, a first two-wire line, one or more four-wire carrier systems transmitted via an orbiting satellite, 
and a second two-wire line, which echoes are reflected back into said receiver by said first and second 
40 two-wire lines respectively, said method comprising the steps of: 

(a) repeatedly transmitting a cyclic signal that has a period that is equal to the duration of said sets 
of encoded frequency-domain samples; 

(b) transforming into the frequency domain those first sets of samples of the signals which are 
received via said near echo path before said far echo has arrived at the receiver of said modem; 

45 (c) dividing the transform of said first sets of samples of step (b) by the transform of said cyclic 

signal to form a first estimate of the frequency-domain characteristics of said near echo path; 

(d) calculating the time-domain parameters of said near echo path by transforming said first estimate 
of step (c) into the time domain, and setting to zero all samples outside the span of said near echo; 

(e) calculating a smoothed set of frequency-domain parameters of said near echo path by 
so transforming said time-domain parameters of step (d) into the frequency domain; 

(f) averaging all sets of said time-domain and frequency-domain parameters of steps (d) and (e) 
respectively that are received before said far echo arrives at said receiver; 

(g) calculating the time-domain samples of the estimated near echo using said averaged sets of 
parameters of step <f); 

55 (h) generating a first set of samples of said far echo by subtracting said samples of step (g) from 

those samples of said actual received signal that are received after said far echo arrives at said 
receiver; 

(I) generating a second set of time-domain samples of said far echo signal by rotating said first set 
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of samples of step (h) by the angles needed to compensate for the estimated frequency offset of 
said far echo path; 

(j) transforming said second set of samples of step (i) into the frequency domain; and 
<k) for all sets of samples of said received signal that comprise both near and far echoes after the 
s first, averaging the change of phase angle of said frequency-domain samples and thereby calculat- 

ing an update of said estimated frequency offset. 

11. A method of initially estimating the parameters of near and far . echoes of a transmitted signal in a data 
modem comprising a multicarrier transmitter which transforms blocks of encoded frequency-domain 

10 samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 
operation, said transmitter and receiver being both connected to a communication link that comprises, 
in turn, a first two-wire line, one or more four-wire carrier systems transmitted via an orbiting satellite, 
and a second two-wire line, which echoes are reflected back into said receiver by said first and second 
two-wire lines respectively, said method comprising the steps of: 

f5 (a) repeatedly transmitting a cyclic signal that has a period that is equal to the duration of said sets 

of encoded frequency-domain samples; 

(b) transforming into the frequency domain those first sets of samples of the signals which are 
received via said near echo path before said far echo has arrived at the receiver of said modem; 

(c) dividing the transform of said first sets of samples of step (b) by the transform of said cyclic 
20 signal to form a first estimate of the frequency-domain characteristics of said near echo path; 

(d) calculating the time-domain parameters of said near echo path by transforming said first estimate 
of step (c) into the time domain, and setting to zero all samples outside the span of said near echo; 

(e) calculating a smoothed set of frequency-domain parameters of said near echo path by 
transforming said time-domain parameters of step (d) into the frequency domain; 

25 (f) averaging all sets of said time-domain and frequency-domain parameters of steps (d) and (e) 

repectively that are received before said far echo arrives at said recever; 

(g) calculating the frequency-domain paramters of the combined near and far echo paths by 
transforming into the frequency domain sets of those samples of said actual received signal that are 
received after said far echo arrives at said receiver; 
30 (h) generating a first estimate of the frequency-domain paramters of said far echo path by 

subtracting from said frequency-domain parameters of said combined echo paths of step (g) said 
smoothed and averaged frequency-domain parameters of said near echo path of step (f); 
(i) transforming said first estimate of step(h) into the time domain; 

(j) generating a second set of time-domain parameters of said far echo signal by rotating said first 
35 set of parameters of step (i) by the angles needed to compensate for the estimated frequency offset 

of said far echo path; 

(k) transforming said second set of parameters of step (h) into the frequency domain; and 
(I) for all sets of samples of said received signal that comprise both near and far echoes after the 
first, averaging the change of phase angle of said frequency-domain parameters and thereby 
40 calculating an update of said estimated frequency offset. 

12. Apparatus for cancelling the echoes of a signal transmitted via a two-wire line in a data modem 
comprising a multicarrier transmitter which transforms blocks of encoded frequency-domain samples 
into blocks of time-domain samples and a multicarrier receiver which performs the inverse operation, 

45 said transmitter and receiver being both connected to said two-wire line, which echoes are reflected 
back into said receiver by an imperfect match between the output impedance of said modem and the 
input impedance of said two-wire line, said apparatus comprising: 

means for estimating the frequency-domain characteristics of the path traversed by said echo; 

50 

means for transforming said estimated frequency-domain characteristics into the time-domain to 
form an estimate of the impulse response of said echo path; 

means for calculating the expected echo signal of said transmitted signal based upon said 
55 estimated impulse response; 

means for subtracting said expected echo signal from the actual signal received from said two-wire 

line; 
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means for operating upon said actual received signal to form a modified received signal of which 
the echo component appears to be cyclic; 

5 means for transforming the time samples of said modified received signal into the frequency 

domain; 

means for correlating said transformed modified received signal with said encoded frequency- 
domain samples of said transmitted signal; and 

70 

means for changing said estimated frequency characteristics of said echo path in order to reduce 
said correlation. 

13. The apparatus of claim 12 further comprising means for performing said transformations from time 
75 doamin to frequency domain using a discrete Fourier transform, and said transformations from 

frequency domain into time domain using an inverse discrete Fourier transform. 

14. The apparatus of claim 12 further comprising means for performing said transformations from time 
domain to frequency domain using a fast Fourier transform, and said transformations from frequency 

20 domain into time domain using an inverse fast Fourier transform. 

15. The apparatus of claim 12, further comprising means for augmenting said encoded frequency-domain 
samples by some other samples that are not used to convey data information, so that the composite 
set occupies a frequency band that is wider than that occupied by said encoded frequency-domain 

25 samples. 

16. The apparatus of ome of claims 12 to 15 further comprising means for smoothing said frequency- 
domain parameters, said means comprising: 

30 means for changing said frequency-domain characteristics in order to generate a first set of new 

characteristics; 

means for transforming said first set of new characteristics into the time-domain, and truncating the 
time-domain impulse response in order to generate a new set of parameters of said estimated impulse 
35 response; and 

means for transforming said new, truncated estimated impulse response back into the frequency 
domain to generate a second set of new frequency-domain characteristics. 

40 17. The apparatus of claim 16, further comprising means for performing said transformations using Inverse 
Fast Fourier Transforms and Fast Fourier Transforms. 

18. Apparatus for cancelling a far echo of a transmitted signal in a data modem comprising a multicarrier 
transmitter which transforms blocks of encoded frequency-domain samples into blocks of time-domain 
45 samples and a multicarrier receiver which performs the inverse operation, said transmitter and receiver 
being both connected to a communication link that comprises, in turn, a first two-wire line, one or more 
four-wire terrestrial carrier systems, and a second two-wire line, which far echo is reflected back to said 
receiver from said second two-wire line, and which may be offset in frequency, said apparatus 
comprising: 

so 

means for forming first estimates of the frequency-domain characteristics and the frequency offset 
of the path traversed by said far echo; 

means for transforming said estimated frequency-domain characteristics into the time-domain to 
55 form an estimate of the impulse response parameters of said far echo path; 

means for calculating the expected time-domain samples of the echo of said transmitted signal 
based upon said estimated Impulse response parameters; 
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means for rotating said expected time-domain samples of by an angle that increased linearly with 
time in proportion to said estimate of the frequency offset; 

5 means for subtracting said rotated expected echo signal of from the actual signal received from 

said two-wire line in order to calculate a modified received signal; 

means for transforming a first block of time-domain samples of said modified received signal into 
the frequency domain; 

w 

means for correlating said transformed modified received signal with said encoded frequency- 
domain samples of said transmitted signal; 

means for calculating a new set of estimated frequency-domain echo path characteristics by 
is decreasing each of said characteristics by a small quantity that is proportional to said correlation; 

means for calculating a set of estimates of the phase angles by which each of said new frequency 
characteristics has rotated; 

20 means for averaging said changes of phase angles over the set of said characteristics; and 

means for filtering said average phase angle to generate an improved estimate of the frequency 
offset of said far echo. 

25 19. Apparatus for initially estimating the parameters of near and far echoes of a transmitted signal in a data 
modem comprising a multicarrier transmitter which transforms blocks of encoded frequency-domain 
samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 
operation, said transmitter and receiver being both connected to a terrestrial communication link that 
comprises, in turn, a first two-wire line, one or more four-wire terrestrial carrier systems, and second 

30 two-wire line, which near and far echoes are reflected back into said receiver by said first and second 
two-wire lines respectively, said apparatus comprising: 

means for repeatedly transmitting a cyclic signal that has a period that is equal to the duration of 
said sets of encoded frequency-domain samples, and longer than the delay of said far echo; 

35 

means for transforming into the frequency domain sets of samples of the signals received via said 
echo paths; 

means for dividing said transform of said sets of samples by the transform of said cyclic signal to 
40 form a first estimate of the combined frequency-domain characteristics of both of said near and far 
echo paths; 

means for calculating the time-domain parameters of said near echo path by transforming said first 
estimate into the time domain, and setting to zero all samples outside the span of said near echo; 

45 

means for calculating a first set of time-domain parameters of said far echo by setting to zero all 
samples of said transformed first estimate which are outside the span of said far echo; 

means for generating a second set of time-domain parameters of said far echo by rotating said first 
so set of time-domain parameters by the angles needed to compensate for the estimated frequency offset 
of said far echo path; 

means for transforming said second set of time-domain parameters into the frequency domain; and 

55 means for averaging, for all sets of samples of said actual received signal after the first, the change 

of phase angle of said frequency-domain parameters and thereby calculating an update of the 
estimated frequency offset. 
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20. Apparatus for initially estimating the parameters of near and far echoes of a transmitted signal in a data 
modem comprising a multicarrier transmitter which transformsblocks of encoded frequency-domain 
samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 
operation, said transmitter and receiver being both connected to a communication link that comprises, 

5 in turn, a first two-wire line, one or more four-wire carrier systems transmitted viy an orbiting satellite, 
and a second two-wire line, which echoes are reflected back into said receiver by said first and second 
. two-wire lines respectively, said apparatus comprising: 

means for repeatedly transmitting a cyclic signal which has a period that is equal to the duration of 
to said sets of encoded frequency-domain samples; 

means for transforming into the frequency domain those first sets of samples of the signals which 
are received via said near echo path before said far echo has arrived at the receiver of said first 
modem; 

;s 

means for dividing the transform of said first sets of samples by the transform of said cyclic signal 
to form a first estimate of the frequency-domain characteristics of said near echo path; 

means for calculating the time-domain parameters of said near echo path by transforming said first 
20 estimate into the time domain, and setting to zero all samples outside the span of said near echo; 

means for calculating a smoothed set of frequency-domain parameters of said near echo path by 
transforming said time-domain parameters into the frequency domain; 

25 means for averaging all sets of said time-domain and frequency-domain parameters respectively 

that are received before said far echo arrives at said receiver; 

means for calculating the time-domain samples of the estimated echo using said averaged sets of 
parameters; 

30 

means for subtracting said time-domain samples from those samples of said actual received signal 
that are received after said far echo arrives at said receiver; 

means for generating a second set of time-domain samples of said far echo signal by rotating said 
35 first set of samples by the angles needed to compensate for the estimated frequency offset of said far 
echo path; 

means for transforming said second set of samples into the frequency domain; and 

40 means for averaging, for all sets of samples of said received signal that comprise both near and far 

echoes after the first, the change of phase angle of said frequency-domain parameters and thereby 
calculating an update of the estimated frequency offset. 

21. Apparatus for initially estimating the parameters of near and far echoes of a transmitted signal in a data 
45 modem comprising a multicarrier transmitter which transforms blocks of encoded frequency-domain 

samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 
operation, said transmitter and receiver being both connected to a communication link that comprises, 
in turn, a first two-wire line, one or more four-wire carrier systems transmitted via an orbiting satellite, 
and a second two-wire line, which echoes are reflected back into said receiver by said first and second 
so two-wire lines repectively, said apparatus comprising: 

means for repeatedly transmitting a cyclic signal which has a period that is equal to the duration of 
said sets of encoded frequency-domain samples; 

55 means for transforming into the frequency domain those first sets of samples of the signals which 

are received via said near echo path before said far echo has arrived at the receiver of said first 
modem; 
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20. Apparatus for initially estimating the parameters of near and far echoes of a transmitted signal in a data 
modem comprising a multicarrier transmitter which transformsblocks of encoded frequency-domain 
samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 
operation, said transmitter and receiver being both connected to a communication link that comprises, 

5 in turn, a first two-wire line, one or more four-wire carrier systems transmitted viy an orbiting satellite, 
and a second two-wire line, which echoes are reflected back into said receiver by said first and second 
two-wire lines respectively, said apparatus comprising: 

means for repeatedly transmitting a cyclic signal which has a period that is equal to the duration of 
70 said sets of encoded frequency-domain samples; 

means for transforming into the frequency domain those first sets of samples of the signals which 
are received via said near echo path before said far echo has arrived at the receiver of said first 
modem; 

means for dividing the transform of said first sets of samples by the transform of said cyclic signal 
to form a first estimate of the frequency-domain characteristics of said near echo path; 

means for calculating the time-domain parameters of said near echo path by transforming said first 
20 estimate into the time domain, and setting to zero all samples outside the span of said near echo; 

means for calculating a smoothed set of frequency-domain parameters of said near echo path by 
transforming said time-domain parameters into the frequency domain; 

25 means for averaging all sets of said time-domain and frequency-domain parameters respectively 

that are received before said far echo arrives at said receiver; 

means for calculating the time-domain samples of the estimated echo using said averaged sets of 
parameters; 

30 

means for subtracting said time-domain samples from those samples of said actual received signal 
that are received after said far echo arrives at said receiver; 

means for generating a second set of time-domain samples of said far echo signal by rotating said 
35 first set of samples by the angles needed to compensate for the estimated frequency offset of said far 
echo path; 

means for transforming said second set of samples into the frequency domain; and 

40 means for averaging, for all sets of samples of said received signal that comprise both near and far 

echoes after the first, the change of phase angle of said frequency-domain parameters and thereby 
calculating an update of the estimated frequency offset. 

21. Apparatus for initially estimating the parameters of near and far echoes of a transmitted signal in a data 
45 modem comprising a multicarrier transmitter which transforms blocks of encoded frequency-domain 

samples into blocks of time-domain samples and a multicarrier receiver which performs the inverse 
operation, said transmitter and receiver being both connected to a communication link that comprises, 
in turn, a first two-wire line, one or more four-wire carrier systems transmitted via an orbiting satellite, 
and a second two-wire line, which echoes are reflected back into said receiver by said first and second 
so two-wire lines repectively, said apparatus comprising: 

means for repeatedly transmitting a cyclic signal which has a period that is equal to the duration of 
said sets of encoded frequency-domain samples; 

55 means for transforming into the frequency domain those first sets of samples of the signals which 

are received via said near echo path before said far echo has arrived at the receiver of said first 
modem; 
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means for dividing the transform of said first sets of samples by the transform of said cyclic signal 
to form a first estimate of the frequency-domain characteristics of said near echo path; 

means for calculating the time-domain parameters of said near echo path by transforming said first 
estimate into the time domain, and setting to zero all samples outside the span of said near echo; 

means for calculating a smoothed set of frequency-domain parameters of said near echo path by 
transforming said time-domain parameters into the frequency domain; 

means for averaging all sets of said time-domain and frequency-domain parameters respectively 
that are received before said far echo arrives at said receiver, 

means for calculating the frequency-domain parameters of the combined near and far echo paths 
by transforming into the frequency domain sets of those samples of said actual received signal that are 
received after said far echo arrives at said receiver; 

means for generating a first estimate of the frequency-domain parameters of said far echo path by 
subtracting from said frequency-domain parameters of said combined echo paths said smoothed and 
averaged frequency-domain parameters of said near echo path; 

means for transforming said first estimate into the time domain; 

means for generating a second set of time-domain parameters of said far echo signal by rotating 
said first set of parameters by the angles needed to compensate for the estimated frequency offset of 
said far echo path; 

means for transforming said second set of parameters into the frequency domain, and 

means for averaging, for all sets of samples of said received signal that comprise both near and far 
echoes after the first, the change of phase angle of said frequency-domain parameters and thereby 
calculating an update of the estimated frequency offset. 
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